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Abstract 

Variational principles for magnetohydrodynamics were introduced 
by previous authors both in Lagrangian and Eulerian form. In a previ- 
ous work Yahalom & Lynden-Bell introduced a simpler Eulerian vari- 
ational principles from which all the relevant equations of magnetohy- 
drodynamics can be derived. The variational principle was given in 
terms of six independent functions for non-stationary flows and three 
independent functions for stationary flows. This is less then the seven 
variables which appear in the standard equations of magnetohydro- 
dynamics which are the magnetic field B the velocity field v and the 
density p. In this work I will improve on the previous results show- 
ing that non-stationary magnetohydrodynamics should be described 
by four functions . 

1 Introduction 

Variational principles for magnetohydrodynamics were introduced by pre- 
vious authors both in Lagrangian and Eulerian form. Sturrock [1] has dis- 
cussed in his book a Lagrangian variational formalism for magnetohydro- 
dynamics. Vladimirov and Moffatt [2] in a series of papers have discussed 
an Eulerian variational principle for incompressible magnetohydrodynam- 
ics. However, their variational principle contained three more functions in 
addition to the seven variables which appear in the standard equations of 
magnetohydrodynamics which are the magnetic field B the velocity field v 
and the density p. Kats [3] has generalized Moffatt 's work for compress- 
ible non barotropic flows but without reducing the number of functions and 
the computational load. Moreover, Kats have shown that the variables he 
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suggested can be utilized to describe the motion of arbitrary discontinuity 
surfaces [U [5]. Sakurai [6] has introduced a two function Eulerian varia- 
tional principle for force-free magnetohydrodynamics and used it as a basis 
of a numerical scheme, his method is discussed in a book by Sturrock pQ. A 
method of solving the equations for those two variables was introduced by 
Yang, Sturrock & Antiochos [7J. In a recent work Yahalom Sz Lynden-Bell 
[8j [9] have combined the Lagrangian of Sturrock [1] with the Lagrangian of 
Sakurai [6] to obtain an Eulerian Lagrangian principle depending on only 
six functions. The vanishing of the variational derivatives of this Lagrangian 
entail all the equations needed to describe barotropic magnetohydrodynam- 
ics without any additional constraints. The equations obtained resemble 
the equations of Frenkel, Levich & Stilman [TO] (see also [H]). Further- 
more, it was shown that for stationary flows three functions will suffice in 
order to describe a Lagrangian principle for barotropic magnetohydrody- 
namics. The non-singlevaluedness of the functions appearing in the reduced 
representation of barotropic magnetohydrodynamics was discussed in par- 
ticular with connection to the topological invariants of magnetic and cross 
helicities. It was shown how the conservation of cross helicity can be easily 
generated using the Noether theorem and the variables introduced in that 
paper. In the current paper I improve on the previous results and show that 
four functions are enough to describe a general non stationary barotropic 
magnetohydrodynamics, the idea is borrowed from [12] see also |13| . I14| , [T5] . 

The plan of this paper is as follows: First I introduce the standard nota- 
tions and equations of barotropic magnetohydrodynamics. Next I introduce 
the potential representation of the magnetic field B and the velocity field v. 
This is followed by a review of the Eulerian variational principle developed 
by Yahalom & Lynden-Bell [8j [9] . After those introductory sections I will 
present the four function Eulerian variational principles for non-stationary 
magnetohydrodynamics. 

2 The standard formulation of barotropic magne- 
tohydrodynamics 

The standard set of equations solved for barotropic magnetohydrodynamics 
are given below: 




dB 




V • B = 0, 
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^ + v-(H = o, (3) 

dv ,dv , . (V x B) x B 
P dt = P{ dt + ( ^ ' V) } = " Vp(/0) + 47 ' (4) 

The following notations are utilized: is the temporal derivative, ^ is 
the temporal material derivative and V has its standard meaning in vector 
calculus. B is the magnetic field vector, v is the velocity field vector and p 
is the fluid density. Finally p(p) is the pressure which we assume depends 
on the density alone (barotropic case). The justification for those equations 
and the conditions under which they apply can be found in standard books 
on magnetohydro dynamics (see for example PQ). Equation ([I]) describes the 
fact that the magnetic field lines are moving with the fluid elements (" frozen" 
magnetic field lines), equation ([2]) describes the fact that the magnetic field 
is solenoidal, equation (|3|) describes the conservation of mass and equation 
@ is the Euler equation for a fluid in which both pressure and Lorentz 
magnetic forces apply. The term: 

is the electric current density which is not connected to any mass flow. 
The number of independent variables for which one needs to solve is seven 
(v, B, p) and the number of equations (jl|3|4p is also seven. Notice that 
equation ([2]) is a condition on the initial B field and is satisfied automatically 
for any other time due to equation ([1]). Also notice that p(p) is not a variable 
rather it is a given function of p. 

3 Potential representation of vector quantities of 
magnetohydrodynamics 

It was shown in [8] that B and v can be represented in terms of five scalar 
functions a, (3, x, Vi v - Following Sakurai [6] the magnetic field takes the 
form: 

B = Vx x Vry. (6) 

Hence B satisfies automatically equation ([2]) and is orthogonal to both Vx 
and Vr/. A similar representation was suggested by Dungey [16] but not in 
the context of variational analysis. The above expression can also describe 
a magnetic field with non-zero magnetic helicity as was demonstrated in [8] . 
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Moreover, the velocity v can be represented in the following form: 

v = W + aVx + PVrj. (7) 

this representation is a generalization of the Clebsch representation |17] for 
magnetohydrodynamics. 

4 The Action of Barotropic Magnetohydrodynam- 
ics 

It was shown in [8j that the action of barotropic magnetohydrodynamics 
takes the form: 

A = f Cd 3 xdt, 



C = -p 

1 ,^ 



Svr (V X x Vr?) 2 , (8) 

in which e(p) is the specific internal energy. Taking the variational deriva- 
tives to zero for arbitrary variations leads to the following set of equations: 

% + v ■ 0*0 = 0, (9) 

l = °- < 10 > 

_ = _ u; (12 
dt 2 v ; 

in which u> is the specific enthalpy. 



p 



(13) 



(14) 



In all the above equations B is given by equation ([6]) and v is given by 
equation ([7]) . The mass conservation equation ([3]) is readily obtained. Now 
one needs to show that also equation ([T]) and equation (jj]) are satisfied. 
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It can be easily shown that provided that B is in the form given in equa- 
tion ©, and equation (fTU|) and equation (fTT|) are satisfied, then equations 
(PQ) are satisfied. 

We shall now show that a velocity field given by equation (|7|), such that 
the equations for a, 3, x, T), v satisfy the corresponding equations (|9I1UI11I121 
I13|14p must satisfy Euler's equations. Let us calculate the material deriva- 
tive of v: 

dv dVv da-. dVx . dVrj 
~dt = ~df ~dt X a ~dF ~dt V (15) 
It can be easily shown that: 

dVv ^,dv ^* dv ^* , 1 _q \ ^, dv 



dt ^ dt ^ Vk dxk 2^ W ^ dx^ 

dt dt dxk dxk 

dVx = ^dX _ y Vk ^]<_ = _y Vk ^2L n6) 
dt dt dx k dxk 

In which x k is a Cartesian coordinate and a summation convention is as- 
sumed. Equations (|10lllll2p where used in the above derivation. Inserting 
the result from equations (|13|14|16j) into equation (fT5l) yields: 

— = V- fc ( — | g dx | 3 dv ) | V( l y 2 w ) 
dt dxk dxk dx k 2 

+ -((Vr?- j)Vx-(V X - J)Vv) 
P 

= -Vv k v k + V{\v 2 -w) + -Jx (Vx x Vry) 
2 p 

= -Yl + ^JxB. (17) 
P P 

In which we have used both equation ([6]) and equation ([7]) in the above 
derivation. This of course proves that the barotropic Euler equations can be 
derived from the action given in equation ([8]) and hence all the equations 
of barotropic magnetohydrodynamics can be derived from the above action 
without restricting the variations in any way except on the relevant bound- 
aries and cuts. The reader should take into account that the topology of 
the magnetohydrodynamic flow is conserved, hence cuts must be introduced 
into the calculation as initial conditions. 
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5 A Simpler Action for Barotropic Magnetohy- 
drodynamics 

Can we obtain a further reduction of barotropic magnetohydro dynamics? 
Can we formulate magnetohydro dynamics with less than the six functions 
a, (3, x, v i The answer is yes, in fact four functions x> rj, u, p will suffice. 
To see this we may write the two equations (11011 1|) as equations for a, {3 
that is: 

^ = ^ + v-Vx = ^- + (Vv + aV X + f3Vr ] ).Vx = 0, 
^ = ^ + # ■ = ^ + (W + aV X + /3Vr/) • Vr? = 0, (18) 
in which we have used equation ([7]). Solving for a, (3 we obtain: 

(V??) 2 (^ + W ■ V X ) - (Vry • V X )(g + W • Vr?) 



(V^-Vx) 2 -^) 2 ^) 
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(V^-Vx) 2 -^^) 2 

Hence a and [3 are not free variables any more, but depend on x,Vi u - More- 
over, the velocity v now depends on the same three variables X , Vi u '- 

v = Vv + a[x,r],v)Vx + f3[x,i], v\Vr). (20) 

Since v is given now by equation (|20|) it follows that the two equations 
(jlOllip are satisfied identically and need not be derived from a variational 
principle. The above equation can be somewhat simplified resulting in: 

B 2 dt dt 1 J 

= fK^-|vr))xB + B(W-B)] (21) 

Hence the velocity v is partitioned naturally into two components one which 
is parallel to the magnetic field and another one which is perpendicular to 
it: 

V = V± + Vu 

v± = -i(^Vx-?rV?7) xB, v„ = M-(Vv-B). (22) 
B 2 dt dt ' 11 B 2 
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Inserting the velocity representation (|2~T|) into equation (fl~9|) will lead to the 
result: 



a 



P 



Vr)-(Bx (v-Vu)) 

Vx-0 x (v-Vv)) 
B~ 2 



(23) 



Finally equations (|19|) should be substituted into equation ([8]) to obtain a 
Lagrangian density C in terms of x, V> v i P- 



C[x,ri,v,P\ 



dv dx dfi 

~ ^(VxxVr?) 2 . 



Using equations (|23p this can be written as: 



£[ Xj v , u, p] = p[-v 2 - — - e(p)\ - —B 2 
" ' FJ Pl 2 dt 8?r 



(24) 



(25) 



were v is given by equation (|21|) and by equation ([U]). Or more explicitly 
as: 



£[x,v,», p] 



1 



2 (V.xV,)^ "t - ^ X Tt + X ^ X ^ 



1^ , (Vx X V??) 2 



(26) 



This Lagrangian density admits an infinite symmetry group of transforma- 
tions of the form: 

f) = f)(x,v), X = X(X,V), (27) 

provided that the absolute value of the Jacobian of these transformation is 
unity: 

9(f), x) 



1. 



d(v,x) 

In particular the Lagrangian density admits an exchange symmetry: 



(28) 



v = x, x = v- 



(29) 
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As a consequence of the double infinite symmetry group we have two local 
conservation laws given by the two equations (jlOllip . Taking the variational 
derivatives of the action defined using equation (|26p to zero for arbitrary 
variations leads to the following set of equations: 

t? + v-(H = o, ( 3 °) 



dt 



-v 2 - w, (31) 



dt 2 

da[x,V,v] = ■ J 
dt p 



(32) 



(33) 



dt p 

Those equations should be solved for x,ry, u, p. Equations (|32|33j) contain 
a complicated linear combination of the second derivatives and ^Ly. 
This is inconvenient numerically therefore the following approach is recom- 
mended. Taking the partial temporal derivative of the two equations (|10lll|) 
we obtain: 

Using the expression ^ from equation (|17p we obtain an explicit expression 
for the second derivatives of the form: 

d 2 y -> -> 1 -» - -> -> dy 

_f = tf-V )v + Vw--JxB ).V X ~v-V^ 
ot z p at 

d 2 n -» -» 1 -» -» -> -> <9?7 

— f = ((tf. V)£ + Vic - -J x B) • Vrj-v- V— -. (35) 

at* 3 p ot 

Hence we arrived at a four function formalism for barotropic magnetohy- 
drodynamics which can be derived from a Lagrangian. Notice, however, 
that this formalism contains two first order equations and two second order 
equations, while our previous six function formalism [8] contained six first 
order equations. 



6 Conclusion 

We have shown that barotropic magnetohydrodynamics can be represented 
in terms of four scalar functions %i V) v t P instead of the seven quantities 
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which are the magnetic field B the velocity field v and the density p. Antic- 
ipated applications include stability analysis and the description of numer- 
ical schemes using the described variational principles, exceed the scope of 
this paper. 

It was shown by the author [T2] that variational principles can be used 
directly for numerical analysis (simulation) without the need to refer to the 
field equations. This mathematical construction may lead to better algo- 
rithms for simulating magnetohydro dynamics in terms of the needed com- 
puter memory and CPU time. This approach was applied to potential flows 
in a series of papers [19] [20]. I21j . Moreover, it was implemented in a user 
friendly software package FLUIDEX (which can be down loaded from the 
web site www.fluidex-cfd.com). A variational formalism of magnetohydro- 
dynamics should serve the same use. 

As for stability analysis I suspect that for achieving this we will need to 
add additional constants of motion constraints to the action as was done by 
[22] . hopefully this will be discussed in a future paper. 
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